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ABSTRACT 



Holmberg IX X-l is an archetypal ultraluminous X-ray source (ULX). Here 
we study the properties of the optical counterpart and of its stellar environ- 
ment using optical data from SUBARU/Faint Object Camera and Spectrograph, 
GEMINI/GMOS-N and Hubble Space Telescope (HST)/Advanced Camera for 
Surveys, as well as simultaneous Chandra X-ray data. The V ~ 22.6 spectro- 
scopically identified optical counterpart is part of a loose cluster with an age 
< 20 Myr. Consequently, the mass upper limit on individual stars in the asso- 
ciation is about 20 M . The counterpart is more luminous than the other stars 
of the association, suggesting a non-negligible optical contribution from the ac- 
cretion disk. An observed UV excess also points to non-stellar light similar to 
X-ray active low-mass X-ray binaries. A broad He IIA4686 emission line identi- 
fied in the optical spectrum of the ULX further suggests optical light from X-ray 
reprocessing in the accretion disk. Using stellar evolutionary tracks, we have 
constrained the mass of the counterpart to be > 10 M@, even if the accretion 
disk contributes significantly to the optical luminosity. Comparison of the pho- 
tometric properties of the counterpart with binary models show that the donor 
may be more massive, > 25M , with the ULX system likely undergoing case AB 
mass transfer. Finally, the counterpart exhibits photometric variability of 0.14 
mag between two HST observations separated by 50 days which could be due to 
ellipsoidal variations and/or disk reprocessing of variable X-ray emission. 

Subject headings: accretion, accretion disks - black hole physics - galaxies: individual: 
Holmberg IX - galaxies: star clusters - X-rays: binaries - X-rays: individuals: 
Holmberg IX X-l 
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Introduction 



Ultraluminous X-ray sources (ULXs) are compact X-ray sources that are located away 
from the nucleus of their host galaxy, emitting well above the Eddington limit of a 20 M Q 
black hole (Lx ~ 3 x 10 39 erg s _1 ) assuming isotropic emission. It is now agreed that 
ULXs are accreting systems powered by black holes, due to the high X-ray variability they 
sometimes display on timescales as short as minutes. Among the remaining questions is 
the nature of these black holes: are they stellar-mass black holes similar to Galactic X-ray 
ligher accretion rates or do they contain intermediate-mass black holes? 



bin aries but with 

See iRoberts! (120071 ) for a recent X-ray review. 



Holmberg IX X-l (HoIX X-l also known as M81 X-9) is a prototype ULX located 



in Holmberg IX, an irregular dwarf galaxy cataloged by 



Holmberg et al. 



(11974 ) belonging 



to the M81 group of galaxies. It is located at about 3.6 Mpc (the distance of M81; 



Freedman et al. 



19941 ). corresponding to a distance modulus m — M = 27.78. Holmberg IX 



may be a tidal dwarf galaxy resulting from the major interaction b etween M81 , M82 



and NGC 3077 that occurred 200 Myr ago or an old dwarf galaxy (jSabbi et al. 



2008|). 



Indeed, this group of gala xies is located in a giant H I cloud with tidal bridges connecting 



the different components (jYun et al. 



19941 ) . Numerous optical knots have been detected 
along the H I filaments where there are overdensities, Holmberg IX being one of the most 
prominent such objects. 



Discovere d by Einstein (IFabbiano 



satellites (see 



Vierdayanti et al 



19881 ) and studied with all the following X-ray 



(120101 ) for an X-ray review), HoIX X-l is an extremely 
luminous X-ray source with an isotropic X-ray luminosity (0.3-10 keV) usually above 
10 40 erg s -1 and thus is one of the most luminous ULXs. The X-ray spectrum of 
this ULX has been interpreted in two different ways, either as a cool accretion disk 
(fcTi n ~ 0.25 keV) with an additional power law or as a slim disk characterized by a hot 
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thermal component (kTi n ~ 1.4— 1.8 keV). This has led to two opposite conclusions 
with the former interpretation suggesting the presenc e of an intermedia te-mass black hole 



(M ~ 1000 M ) accreting at a sub-Eddington rate (jMiller et al. 



200J) and the latter 



i mplying a stel 



(ITsunoda et al 



ar-ma ss black hole (M ~ 10 M Q ) accreting at a super-Eddington rate 



2006). Recent 



corona ([Gladstone et al. 



2009 



y, a model consisting o f a disk plus a cool, optically thick 



Vierdayanti et al. 



20101 ) has been applied to high-quality 



X-ray spectra of ULXs and in particular to HoIX X-l. This model suggests a stellar-mass 
black hole with super-Eddington accretion. 

HoIX X-l is located about 2' (~ 2.1 kpc) northeast of the center of Holmberg IX and 
about 1' (~ 1 kpc) outside of what we can consider the limit of the stellar envelope of the 
galaxy. This ULX is thus present in a zone wher e there is no intense s tellar formation, but 



we note the existence of a close molecular cloud flBrouillet et al. 



19921 ) at 4 0" southeast o f 



its po sition. The Galactic interstellar extinction is low, E(B — V) ~ 0.08 ( ISchlegel et al. 



19981 ). making this source a very good target for an optical photometric and spectroscopic 
study. 

The striking charact eristic of HoI X X-l at optical wavelengths is the presence of a 



huge ionized nebula, see 



Miller! (119951 ) and Figure [D with a size of 300 x 470 pc at the 



distance of 3.6 Mpc. This and further work (jPakull Sz Mirioni 



2002 



20031 ) has shown that 



shocks should play a key role in the nebula even if photoionization from the X-ray source 



is po ssibly required to understand all the properties of the nebula flAbolmasov fc Moiseev 



20081 ) . The optical counterpart of the ULX is the most luminous object belonging to a small 



stellar association located inside the nebula flGrise et al. 



2006 



Ramsey et al. 



20061 ). This 



paper concentrates on the properties of the optical count 



revising and expanding our previous results (IGrise et al 



erpar t and its stellar environment, 



2006|) 
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2. Observations and data analysis 
2.1. SUBARU 

The SUBARU data set comes from two observing programs (Pis: T. Tsuru, M. 
Pakull) and makes use of the SUBA RU/Faint Object Camera and Spectrograph (FOCAS) 



instrument (IKashikawa et al. 



20021 ). This instrument has been used to study the 
counterparts and environment of several ULXs, one of them being HoIX X-l. This data 
set consists of imaging through broadband filters as B, V, and R and through narrowband 
filters as H Q and [O III], as well as spectroscopy. All the raw frames were reduced under 
ESO MIDAS using standard procedures, i.e., correcting for bias and flat fields. The two 
spectra were combined and then wavelength and flux calibrated. 

The local environment to the ULX is not completely resolved, so point-spread function 
(PSF) fitting photometry is necess ary to separate blended objects. We used the package 



DAOPHOT II flStetson et al. 



19901) under ESO MIDAS and we employed a quadratically 



varying PSF in the field. 

Our absolute photometry in the B, V, and R bands uses three consecutive observations 
made on 2003 Jan 26 (see Tabled] for details of the SUBARU observations). Photometric 
transformations were derived with the only observed standard field, Rul52, for which we 
used both the Landolt standard stars and also a set of stars defined by Peter Stetson []. 
This permitted us to obtain a more homogeneous calibration on the range of photometric 
colors and also to use stars fainter than the 14th mag. 

To check our photometry, we used the Sloan Digital Sky Survey (SDSS) data release 6 



Photometry of this secondary set of stars is available on the CADC Web site, 



http:/ /www3.cadc-ccda.hia-iha.nrc-cnrc.gc.ca/community/STETSON/standards/ 
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( lAdelman-McCarthy et al.ll2008f ) which imaged fields around our target. Thanks to its 2.5 m 
telescope, photometry with an accuracy of 2%-3% is possible down to 21-22nd magnitude. 
The filters' system used by the SDSS is noticeably different from the Johnson- Cousins 
system, so we had to transform the u, g and r magnitudes into B, V, R. For that we used 
the empirical transformations of Luptonil (2005) who made a correlation between SDSS 
photometry and Stetson photometry. A possible dependence on color or metallicity related 
to the spectral type of stars has not been taken into account. As our goal is to make a 
global comparison between the two photometric data sets and not an analysis of SDSS 
data, we can content ourselves with this approach. 

The results of this comparison are shown in Figure [2j We note the good linear 
agreement in the three filters. The scatter at faint magnitudes is consistent with the 
photometric errors of the SDSS and is more evident in the R filter, where the errors are 
higher than in the other filters. However, there is an offset in each filter, with the SUBARU 
magnitudes being fainter than the SDSS values (+0.23 mag in B, +0.30 mag in V, and 
+0.22 mag in R). The reason for these offsets are unclear. In the following, all SUBARU 
measurements have been corrected by these offsets. 

Long-slit spectroscopy has been done with the slit centered on the assumed counterpart. 
The observational settings include a 150 lines mm -1 grism, an SY47 filter, and a slit of 0'/6 
in width. This translates into an accessible wavelength range of 4600-7600 A, a dispersion 
of 2.8 A pixel -1 and a resolution A/AA = 500 at 6500 A, i.e., the spectral resolution is 
approximately 8 A. 



See http://sdss.org/dr6/algorithms/sdssUBVRITransform.html 
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2.2. HST 



The second data set used in this study comes from the Hubble Space Telescope 
(HST)/ Advanced Camera for Surveys (ACS) archive (Program GO-9796 ; PI: J. Miller). 
Observations were made on 2004 February 07 with the Wide-Field Camera (WFC) in the 
FA35W, F555W, and F814W filters, and also with the High Resolution Camera (HRC) in 
the F330W filter. One additional observation was made on 2004 March 25 with the WFC 
in the F555W filter. More details about these observations are given in Table [2J We have 
re-analyzed the data using the latest calibration files available at the time of the analysis 
(STSDAS v3.5) and have appli ed the standard correct ions, especially about the spatial 



distortion (Multidrizzle v2.7.2; 



Koekemoer et al. 



20021 ). As in the SUBARU observations, 



the field of int erest looks moderat ely crowded in the HST/ ACS image so we used again 



DAOPHOT II flStetson et al. 



19901 ) for the photometric analysis. This analysi s is similar 



to our previous study of the ULX NGC 1313 X-2 and details can be found in 



Grise et al. 



( 12008 



G08 hereafter). 



2.3. Comparison between HST and SUBARU photometry 

Optical variability is an important parameter in studies of X-ray binaries, since it can 
constrain the orbital period of the system or quantify a possible variability linked to the 
accretion disk. To study photometric variability using different telescopes, we have to make 
sure that the magnitudes derived from the HST emd SUBARU data sets are consistent. As 
for our previous study (G08), we have isolated bright point-like sources in the two data sets 
and compared their luminosities to study the photometric accuracy of our two data sets. 

Figure [3] shows the magnitude differences between FOCAS and ACS in the B and V 
filters for observed stars. For 50% of them, the differences are consistent with the combined 
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photometric errors of the two instruments. For the other objects, the differences are 
almost always below 0.1 mag. Standard deviations are Ob = 0.06 mag and ay = 0.07 mag 
which turn out to be the same as found in our previous study of NGC 1313 X-2 between 
HST/ ACS and ESO/VLT data (G08). We note that deviations of the order of 0.1 mag are 
easily explained by, for example, a 23 mag object blended with a 26 mag object. 

Photometric study of the ULX counterpart with SUBARU, Table |3l was done on 2003 
Jan 26 with three consecutive observations in the three photometric bands B, V and R, 
strongly limiting the risk of seeing intrinsic variability which could have introduced a bias 
in colors estimation. Comparing with HST data, Table HI we note excellent agreement in 
the (B — V) colors with magnitude differences being explained by variability. 



2.4. GEMINI 

The last optical data set used in this study comes from the GEMINI /GMOS-N 
archive (Program ID GN-2008A-Q-49; PI: J. Gladstone). The observations consist of six 
spectra taken consecut ively in the night of 2008 March 07 with the Gemini Multi-Object 



Spectrograph (GMOS) ( iHook et al. 



20041 1 installed on the 8.1 m Gemini North telescope. 
The B600 grism was used, covering a wavelength range of 3160-6070 A with a dispersion of 
0.90 A pixel -1 including a 2 x 2 binning, and the slit width was 0'/75. 

Data reduction was performed using the latest GEMINI IRAF package vl.10. All 
standard tools to reduce long-slit spectroscopy are included and easily handle the three 
CCDs frames. Flat fields were created and used with the processed bias to correct for 
instrumental defects. After wavelength calibration, the six spectra were co-added using the 
ccdclip rejection algorithm to remove cosmic rays and other bad pixels. Careful extraction 
of the spectrum was done due to the trace which appeared to vary by six pixels from one 
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side of the frame to the other. We did not perform an absolute flux calibration. 



2.5. Chandra data 



Simultaneous X-ray and optical observations are quite rare for ULXs, so we took the 
opportunity that two X-ray pointings were done at the same time as the HST observations. 
We retrieved the data from the Chandra archive which consist of two 5 ks ACIS-S 
observations (PI: J. Miller, see details in Table We used the level 2 event files but were 
faced with severe pile-up due to the combination of the high count rate of the source and 
the frame time used, 1.8 s. St andard extra ction of the data l ed to highly d istorted spectra, 



and use of the pile-up model ([Davis 



200 ll ) in Xspec 12.6.01 flArnaud 



19961 ) was not able 



to solve this problem. Thus, we extracted only the counts coming from the PSF wings 
of the source, where pile-up should be negligible. To do so, we used the recipe devised 
by M. Tsujimotclfl that provides a method to correct the Ancillary Response File when 
using events from the outskirts of the PSF. To maximize the signal-to-noise ratio without 
including a large fraction of piled-up events, we used a series of different extraction annuli, 
then fitted an absorbed power law in Xspec and plotted the power-law index versus the 
radius of the inner extraction annulus. The power-law index is unphysically hard (r ~ 0.5) 
at small inner radii, but increases at larger radii until reaching a plateau. This is where 
pile-up becomes negligible and corresponds to an inner extraction radius of 3.5 and 4.0 
pixels for the first and second observation, respectively. The outer radius for the extraction 
was taken as 20 pixels for both observations. 

Results of the spectral fits, involving 286 and 240 counts, respectively, are shown in 
Table [5J where the spectra were grouped to have bins of at least 20 counts. We were 



3 http:/ /www. astro. isas.jaxa.jp/~tsujimot/arfcorr. html 
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unable to use models more complex than an absorbed power law due to the loss of counts 
resulting from the pile - up. Our fits, using a co lumn density fixed to = 1.0 x 10 21 cm -2 



([Gladstone et al. 



2009 



Vierdayanti et al 



20101 ) are quite similar for the two observations 
with a power-law index of ~ 1.8 and an absorbed flux of ~ 1.1 x 10~ n erg s" 1 cm" 2 . This i s 



consistent with the apparent permanent hard spectrum of Ho IX X-l ( IKaaret Sz Feng 



20091 ). 



3. Results and discussion 



3.1. Identification of the optical counterpart 



3.1.1. Photometric identification 



The HST images were individually astrometrically re-calib rated using 10 stars from 



20061 ). The rms of these fits are 



the Two Micron All Sky Survey catalog (ISkrutskie et al. 
between 0'.'14 and 0'.'18. In the Chandra observations, the ULX is located almost on-axis, 
with an offset of only ~ 30". In these conditions, the Chandra X-ray absolute position has a 
90% uncertainty of 0'.'6. We have compared the X-ray positions of some sources outside the 
field shown in Figure S] with the HST/ ACS images and found that they are coincident with 
extended objects, hence preventing a re-calibration. Taking the optical centroids of these 
objects, we found that two of them are consistent with the Chandra X-ray position within 
0'/15 and the third one shows a larger offset, 0'/35. The positions of the ULX in the three 
HST/ ACS observations are all within 1 pixel, i.e., ff!05 and are totally consistent with the 
Chandra position (Table E]). 

We also used for comparison the XMM-Newton position, Table [61 coming from a long 
observation (PI: Tod Strohmayer, ObsID: 0200980101, ex posure time: 100 ks ) and listed 



in the second XMM-Newton serendipitous source catalog ( jWatson et al. 



20091 ). The XMM 



position has a statistical error of 0'.'08 but is affected by a 0'.'35 systematic uncertainty 
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(ler). Thus, the 90% confidence radius, including the optical astrometric uncertainty, is 0'/8. 
Figure S] shows the Chandra and XMM-Newton positions of HoIX X-l with, respectively, a 
0'/67 and a 0'/80 radius error circle (90% confidence, including the error on the calibration 
of the optical image) on an F435W HST/ ACS image, which does not allow one to identify 
a unique counterpart. Specifically, the brightest object is contaminated by a nearby faint 
star, V ~ 24.5, at ~ Cf!2 northwest of its position. 



3.1.2. Spectroscopic identification 

Our SUBARU spectroscopic observations were centered on the presumed counterpart, 
i.e., the brightest object of the small association of stars. The ULX is embedded in the 
optical nebula MH9/10 but fortunately, the central position of the counterpart is coincident 
with a region where the nebula has low intensity. 

The one-dimensional SUBARU spectrum of the most luminous object in the X-ray 
error circle with the slit oriented east-west is shown in Figure [5j The main feature of the 
spectrum is the presence of the A4686 He II emission line, which is a typical signature of 
X-ray binaries and thus identifies this object as the optical counterpart of the ULX. The 
line is spatially coincident with the continuum of the counterpart, Figure [6l and excludes 
the possibility that the emission comes from the nearby fainter star. The He II line is 
not spectrally resolved in the FOCAS spectra. Its equivalent width (EW) is 12 ± 3 A. 
The emission coinc ident with the continuum of the coun terpart is unresolved, Figure El as 



in NGC 1313 X-2 (IPakull et al. 



2006 



Grise et al. 



20091 ). and thus probably reflects the 



emission of the ULX system. 

The GEMINI spectrum confirms this result, see Figure [7J In this case, the slit was 
oriented north-south so the brightest object and nearby fainter star are located at almost 
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the same position along the slit. The He II line is again spatially coincident with the 
brightest object. Thus, the fainter object cannot be the ULX counterpart, in agreement 
with the SUBARU spectrum. The GMOS spectrum has better resolution (pixel scale of 0.9 
versus 2.8 A pixel -1 and spectral resolution of 3.7 versus 8.0 A for SUBARU spectrum), 
allowing the emission lines to be resolved. In the co-added spectra, the He II line has an 
intrinsic FWHM of 7.0 A or 450 km s _1 compared to the nebular lines like [O III]A5007 
with an intrinsic FWHM ~ 3 .1 A or 200 km s -1 . The width of th e He II line is comparable 



with that of NGC 1313 X-2 flPakull et al. 



2006 



Grise et al. 



20091 ) and indicates that the 



He II emission is produced in the ULX syste m rather than bein g of nebular origin. Contrary 



to the behavior observed in NGC 1313 X-2 (1 Grise et al. 



2009|), the EW of the He II line 



is stable in HoIX X-l between the SUBARU and the GEMINI observations, with a value 
of 11 A in the latter which is consistent with the high and persistent X-ray luminosity of 
the ULX. We note that the EW of the He IIA4686 line is more than 10 times higher in this 
object than in the brightest Galactic (and Magellanic) high-mass X-ray binaries which is 
consistent with the higher X-ray luminosity of ULXs. 

The radial velocity of the He II line is similar in the six individual GEMINI observations 
and does not show significant variation between the SUBARU and the GEMINI observations 
with a velocity consistent with those of the nebular lines. Although this result could just 
be a matter of coincidence, we would expect a radial velocity shift if the He II line is indeed 
emitted in the accretion disk and the black hole is not ve r y mas sive. We note that a recent 
dedicated spectroscopic campaign done by 



Roberts et al. 



(120101 ) has shown that the He II 



line displays non-periodic radial velocity variations with rapid variability of the line profile 
on timescales shorter than a day, casting doubt on the possibility to derive a mass function 
for this object with this method. 



Thus, the spectrum of the counterpart of HoIX X-l is similar to other spectra of ULX 



13 



counterparts (e.g, 



Grise et al 



2009 



Kaaret fc Corbel 



20091 ). specifically a blue continuum 



without any spectral features apart from the He II emission line. The two-dimensional 
spectrum of the counterpart is contaminated by the presence of nebular emission lines. 
As the intensity of the nebula is highly spatially variable, we extracted the spectrum by 
subtracting a nebular background taken from a small region close to the ULX counterpart 
on the faint side of the nebula. We note the presence of emission lines such as [O III], 
[O I], H a , and [S II], which are probably residual emission of the nebula due to its highly 
non-uniform intensity. 



3.2. Stellar association around the ULX 

The multicolor SUBARU BVR image (Figure [T]) shows that the environment around 
the X-ray source is dominated by a compact group of blue stars located inside the H Q 
nebula, the R filter also includes the H Q line. Outside of the nebula, we observe a mix of 
blue and red stars, and also foreground stars. On a large-scale HST image (Figured]), we 
note that HoIX X-l is located about 1' east of an area rich in blue stars and is outside the 
standard D 2 5 isophotal radius of the galaxy Holmberg IX. On the contrary, we observe few 
stars north of the ULX, with a diffuse background of red stars close to the limit of detection 
and a few blue stars. 

The ULX clearly does not belong to an intense star-forming region but rather appears 
outside the main body of the galaxy. The tight stellar association around the ULX is 
resolved into some seven stars (Figure [T|) in the SUBARU data thanks to the excellent 
seeing during the observations (Table [1]). The multicolor HST image (Figure |9]) resolves the 
stellar association close to the ULX into dozens of members. 
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3.3. Extinction 



The quality of the photometry and the good agreement between the two data sets are 
two major points. Another key parameter for the interpretation of the data is the extinction 
alon g our line of sight to ward HoIX X-l. The Galactic reddening is low {E(B — V) — 



0.08; 



Schlegel et al. 



19981 ). and studies of the stellar po pulation in Holmberg 



Makarova et al. 



IX have never 



20021 ). To check the 



clearly identified any substantial, local reddening (e.g., 
extinction, we plotted the spectral energy distributions (SEDs) of all the stars located near 
the ULX dereddened by the Galactic extinction, see Figure [10J Given the magnitudes and 
colors of most of those stars (B — V ~ 0.1, V — I ~ 0.0, —5 < My < —2), they have to be 
main-sequence early-type stars, such as 08 to B2 stars, although we cannot rule out that 
a few may be evolved. Comparison of these st ars with SED templates of early-type stars 



( iSchmidt-Kaler 



1982 



Straizys 



1992 



Cox 



2000( 1 show that their slope is not in agreement. 



The extinction must be increased to E(B—V) ~ 0.3 mag to make the observed fluxes 
consistent with the SED templates, Figure EH We note that this estimate is made using 
only the BVI photometric bands (and close related HST bands FA35W , F555W , F814W) 
because the U band is different from the F330W band, the former lying on top of the Balmer 
jump and the latter being located shortward of the jump. We also studied the variation of 
the Balmer decrement (H Q /TLj) in the spectrum of the nebula and we found that MH9, the 
east side of the nebula (along the slit) shows a higher value compared to MH10 (the west 
part), i.e., E(B-V) = 0.26 ± 0.04 mag compared to E(B-V) = 0.09 ± 0.05 mag. This 
suggests differential extinction in the nebula. In addition, the SEDs of several stars located 
in a small region west of the ULX indicate higher reddening, up to E(B — V) ~ 0.8 mag. In 
the remainder of this paper, we will use our best estimate of the extinction for the ULX, 
E(B — V) = 0.26 ± 0.04 mag, with the caveat that there may be some variation of the 
extinction. 
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3.4. Interpretation of the Color-Magnitude diagrams 



Color- magnitude diagrams in B — V and V — I (Figure [TTj) show the predominance of a 
blue main sequence (B — V ~ —0.2, V — I ~ —0.2) but also the existenc e of a g iant branch, 
especially obvious in the V — I diagram, as pointed out by ISabbi et al.l (120081 ) with their 
observations centered on Holmberg IX. The population of stars apparently associated with 
the ULX is nicely resolved with HST and we can see 10 members well identified by their 
positions internal to the nebula. The UV emission as seen in the F330W filter, see Figure HJ 
provides a further clue to their relatively young age. In our study of NGC 1313 X-2 (G08), 
comparison with the F330W image allowed us to discriminate between young stars and 
old field stars in a semi-quantitative manner. Here, the situation is more complex since 
the dominant population is young, likely because the encounter between M81, M82, and 
NGC 3077 triggered a wave of star formation within the past 200 Myr. Nevertheless, we 
can assign the status of "candidates of the association" to all stars visible in the UV image 
and which are visually in the dense association. 

This group of stars does not differ from the other field stars but fits very well in the main 
sequence defined by the other stars. Figure [IT] shows a color-magnitude diagram corrected 
by the assumed redd ening E(B — V) = 0.26 an d an extinction bas e d on the Cardelli la w 



(ICardelli et al. 



19891 ). The Padova isochrones (IMarigo et al 



used are for a metallicity Z = 0. 008, as derived 



of an H II region (region 8 from 



3y 



Miller &: Hodge 



2008 



Girardi et al. 



Makarova et al 



2010|) 



(120021 ) by spectroscopy 



1994; located 2.5' from HoIX X-l) and 



consistent with the dwarf nature of the galaxy. Concerning the age determination, we are 
facing a situation where most of the stars are located near the 1-20 Myr isochrones on 
the main sequence. The width of the main sequence, explained by photometric errors and 
probably some differential reddening, is quite large compared to the separation of the 1-20 
Myr isochrones. We note that there are three stars (five if we consider the probable two 
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reddened stars) with an absolute V magnitude between -4.5 and -5.0. These stars may be 
late-0 type stars or early-B stars evolving toward the supergiant phase. Thus, the age of 
this stellar association is uncertain but is most likely younger than 20 Myr. 

We note that there is a population of stars in the field that is probably older, ~ 50 Myr 
old, as seen by the ~ 10 red supergiants at {My ~ —5; B — V ~ 1.4, V — I ~ 1.3). There 
are no such red supergiants associated with the stellar association around the ULX, but as 
we remarked in G08, the ratio of blue to red supergiants observed in metal-poor galaxies 



(ILanger fc Maeder 



19951 ) being ~ 3 would imply only one or two red supergiants in the 



association which, given the statistics, is not in contradiction. 



A previous study done by iRamsey et al.l (120061 ) found an age of 4-6 Myr using the 
same HST data in the B and V filters. They in fact assumed that the blue stars in the 
color-magnitude diagram are reddened (E(B — V) = 0.2) main-sequence O stars. Even 
though this assumption had no real basis, it seems that our finding, based on a more careful 
extinction estimate, may be consistent with a close value. 

In terms of mass, the Padova evolutionary tracks show that the most luminous stars 
are compatible with a zero-age mass of the order of 10-20 M (Figure fT2l) . We note also 
that in principle, the F330W band would be more discriminating in terms of age and mass 
of the stellar population, but unfortunately there is no available observat ional calibration 



20051 ). Finally, using a 



between the F330W and the more common U band flSirianni et al. 
solar metallicity [Z = 0.019) instead of the sub-solar one assumed here (Z = 0.008) does 
not influence our results in any significant manner. 
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3.5. Properties of the stellar association 



The association of stars around HoIX X-l looks more like an OB association than 
a cluster linked by gravitation. This is something we also saw in NGC 1313 X-2 (G08) 
The integrated magnitudes of the stellar association are Mb ~ —7.6, M y ~ —7.4 and 
Mjjy —7.3. Comparing these magnitudes with Starburst99 simulations (ILeitherer et al. 



19991 ). and assuming that the age of the cluster is 10/20 Myr with a Salpeter initial mass 
function (IMF) or an "a = 3.30" IMF (with masses between 1 and 100 M ), the mass of 
the cluster is about M = (1.5 ± 0.3)/(2.7± 0.3) x 10 3 M Q , which is similar to the association 
surrounding the ULX of our previous study (G08). This group of stars is clearly typical for 
OB associations of the Local Group, both in terms of mass and of size, with its diameter of 
90 pc (~ 5"). This isolated star formation is probably linked to the tidal forces between M81 
and M82 which have created many stars concentrations, looking like eith er dwarf galaxie s 



(like 



lolm berg IX) or less massive associations as recently discovered by 



de Mello et al. 



f l2008h and 



Davidge! fj2008h . 



4. Nature of the ULX optical counterpart 

4.1. Constraints from photometry 

We have shown that the stellar population around the ULX has an age of < 20 Myr 
and an individual upper mass limit of 15lg M & . In contrast to the NGC 1313 X-2 (G08), 
the ULX of our present study has a counterpart brighter than the other stars belonging 
to the stellar association, by almost 1 mag. If the counterpart is a normal star (as 
opposed to the case where the accretion disk may play a predominant role in the optical 
emission), its magnitude My ~ —6.0 would be consistent with an 03-04 star on the main 
sequence, a late-0 type evolved star, or a late-O/early-B lb supergiant. Its dereddened 
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color (B — V)o ~ —0.25 (Table [7j) would be consistent with most of those scenarii, but 
its (V — J)o ~ —0.06 is somewhat too red to be compatible with such interpretation. 
Indeed, the counterpart shows a slightly "red" excess (compared to the other stars of 
the association) since its (V — I) color would be more typical of a late-B star on the 
main-sequence, or a mid supergiant B star. But perhaps the most important feature we 
see is a significant excess in the UV (compared to the other stars of the association). This 
can be seen in the color-magnitude diagram (U — B) versus B, Figure [131 where the ULX 
counterpart is clearly bluer than any other star of the stellar association. Unfortunately, 
this diagram is only illustrative of this UV excess since the magnitudes and colors U, B and 
[U — B) derived from HST/ACS data are not compatible with those used in the available 
evolutionary tracks and isochrones and thus cannot be used to derive age or mass properties 
of the stars. 

The SED of the ULX counterpart, Figure HH leads to the same conclusion. Compared 
to other stars of the association, including the brightest ones in the UV, we can definitely 
see that the ULX counterpart shows a clear excess in the UV and a slight excess in the 
red part of the spectrum. In fact, it shows a 3363-4317 A slope steeper than a 05 star, 
pointing to a very hot continuum in the UV. Figure [14] also shows the SED of the ULX 
NGC 1313 X-2 (G08); NGC 1313 X-2 looks like a normal OB star which is relatively 
unexpected since we have shown in G08 that the optical variability of the source was 
likely associated wit h some X-ray reprocessin g in the system. This deserves more detailed 



study. Interestingly, 



Kaaret fc Corbel! ( 120091 ) showed that the ULX NGC 5408 X-l also 



displays a UV excess at wavelengths shorter than the Balmer edge. They interpret this 
effect as coming from the accretion disk and the photoionized nebula surrounding the 
counterpart. In our case, the main nature of the nebula is quite d i fferent since it seems to 



be mainly shock-ionized (IGrise et al. 



2006 



Abolmasov et al 



2007 



Pakull fc Grise 



2008|) 



but X-ray/UV photoionization is possibly required to explain all the properties of the ULX 
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nebula (lAbolmasov Sz Moiseevl 120081 ). Also, Roberts et al. (private communication) have 
shown, using very deep optical spectra of these two ULXs, that NGC 1313 X-2 may display 
some absorption features coming from the donor star (specifically, a Balmer jump), while 
HoIX X-l shows only a very blue continuum consistent with our findings. In any case, the 
absence of clear absorption lines in both ULXs point to a contribution from the accretion 
disk, which is probably acting at different levels and/or at different wavelengths in both 
sources. 



The X-ray to optical flux ratio, £ = B + 2.5 logF x (e.g., lKaaretll2005l ). with F x the 
2-10 keV unabsorbed X-ray flux density in fiJy and B the d ereddened B magnitude, is 



£ = 2 0.6. This is consistent with the values found in LMXBs (Ivan Paradijs &: McClintock 



19951 ) where the optical emission is dominated by X-rays reprocessing. We remark that if 
the accretion disk and the secondary star had equal luminosities, subtracting the light from 
the accretion disk would correct the position of the counterpart to be consistent with a 
My ~ —5.3 mag star. Thus, it would simply be as luminous as the brightest stars in the 
tight association. The UV excess, as in low-mass X-r ay binaries (LMXBs), would be d ue 



to reprocessed energy radiated in the UV and optical (Ivan Paradijs &: McClintock 



19951). A 



key problem in ULXs is that we do not know the spectral type of the companion star. 



A few models predicting the optical properties o 



literature (IRappaport et al 



Patruno fc Zampieri 



2008 



2005 



Copperwheat et al. 



ULXs haye appeared recently in the 



2005 . 



2007 



Madhusudhan et al 



2008 



Although their conclusions are not unanimous, one common 
result of these studies is that the companion star should be massive enough, M > 10 M , 
to sustain a stable accr etion rate and to allow for X-ray luminosities of 10 40 erg/s. This 



was also pointed out in 



Podsiadlowski et al. 



(120031 ) where the authors show that steady 



accretion can be achieved with a donor star of 8-15 M Q around a 10 M & black hole and 
can lead to X-ray luminosities above 10 40 erg s _1 during 9-5 x 10 5 yr, respectively. These 
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phases of high X-ray luminosity happen when the donor become s a giant and the evolutio n 



is driven by the nuclear evolution of the hydrogen-burning shell ( jPodsiadlowski et al 



Observational resu 


ts also 


2004 




Kaaret et al. 


2006 





Feng fc Kaaret 



companions for some ot 



2008 



Motch et al. 



rer ULXs (e.g., 



2010|)- 



2003). 



Liu et al 



The most recent pub 



( IPatruno fc Zampieri 



2008 



i shed model predicting the optical output of ULXs 



20101 ) is also the only one that takes into account the 
fact that not all the mass lost by the donor star is accreted by the black hole. Instead, 
and especially when the donor star is leaving the main sequence, the rise in the mass 
transfer rate above the Eddington limit implies that there is substantial mass loss from 
the binary that does not contrib ute to the X-ray luminosity and to the reprocessing 
( IPatruno &: Zampierill2008l |2010| ). One of the main differences with previous models is that 
accretion disks cannot be very bright in the optical unless the black hole is massive enough 
and/or the donor star is in the fast expanding giant phase. 

If we apply our photometric results to this model, HoIX X-l is consistent with a 
30-50 M donor around a 20-100 M Q black hole undergoing case A(B) mass transfer. The 
donor is likely to be ascending the giant branch, which would give rise to a larger mass 
accretion rate and thus, a larger X-ray luminosity. However, this sta ge is likely to be very 



short for such a massive star, only a few tens of thousands of years 
20081) . Based on the same arguments applied to NGC 1313 X-2 by 



Patruno fc Zampieri 



Patruno fc Zampieri 



(120 101 ). case B mass transfer is 



a nebula that is ~ 1 Myr old (jPakull fc Grise 



unlikely due to the fact that the ULX is surroun ded by 



2008 



Abolmasov fc Moiseev 



20081 ). Thus, 



case AB, where mass transfe r starts during the main s equence, is more likely since it 



allows several millions years ( IPatruno &: Zampieri 



20081 ) to inject the required energy to 



the nebula. In that framework, our values of V magnitude and B — V color do not allow 
discrimination of the mass of the black hole. Addition of other photometric bands in binary 
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models may help in further constraining the mass of the components. Obtaining a full 
SED of the ULX from the far-UV to the near-infrared would also be helpful in trying to 
disentangle the contributions coming from the accretion disk versus the donor star. 

From our own results, the firm conclusion that we can draw about the companion 
star is that it is likely a quite massive star, M > 10 M , when we take its position in 
the color-magnitude diagrams as those of a non-irradiated star. The same conclusion 



y contamination of the opt ical flux by the accretion 



holds if we take into account a like 
disk. Unfortunately, applied to the iPatruno &: Zampieril (120081 ) model, our results do not 
favor either the "standard" stellar-mass black hole (~ 10-20 M ) or the more massive 
(~ 100 M ) black hole scenario. However, we showed that the total stellar-mass in the 
young stellar association around the ULX is < 10 4 M Q . Thus, HoIX X-l is clearly not 
in a super-star-cluster. This means that whatever the nature of t he compact obje c t, it 



was not formed via a runaway-co a 



and 



escence scenario as proposed by 



Ebisuzaki et al 



(1200 lh 



Portegies Zwart fc McMillan! (120021 ). If the black hole has a mass of the order of 



< 100 Mq, it was probably formed via a more conventional stellar evolutionary scenario. 



This may be possib 



Zampieri fc Roberts 



e tha nks to the sub-solar metallicity of the ULX environment, see 



2009 



for such arguments. 



4.2. Photometric variability 

We have tried to characterize the possible optical variability of the counterpart. We 
have four SUBARU images in the V band, but only two of them have a good enough seeing 
to separate, via PSF fitting, the luminous counterpart from the fainter object which tends 
to contaminate the optical emission. In these two images, the counterpart has a constant V 
magnitude within the photometric errors (22.71 ± 0.04, 22.67 ± 0.04). 
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On the contrary, the two HST/ ACS images in the 1/ band show variability, 
AV = 0.136 ± 0.027 clearly significant at 5er. The amplitude is similar to that of 
NGC 1313 X-2 (G08). This optical variability may originate from ellipsoidal modulation 
reflecting the orbital movement of the donor star, like in massive X-ray binaries. It may 
also come from reprocessing of variable X-ray emission by the accretion disk and X-ray 
heating of the secondary as observed in LMXBs. A combination of both effects is likely, 
since the ULX environment is composed of moderately massive stars with M max ~ 2OM , so 
that the optical emission of the donor is not negligible as in LMXBs, while the ULX X-ray 
output is 10-100 times that of LMXBs, suggesting a non-negligible optical contribution 
from reprocessing. 

The SUBARU data do not permit a study without ambiguity of the intrinsic variability 
of the ULX, but the magnitude V (22.71 ± 0.038) from the 2003 January 26 observation 
is consistent with the range seen by HST. The B — V color is consistent between the two 
observations of the two telescopes (Tables [3] and HJ). We can compare the total magnitude of 
the counterpart and of the nearby object between the two observatories for all observations 
(Table [8]). The agreement is very good between our SUBARU measurements (where we can 
resolve the counterpart from its nearby companion) and the HST observations. The second 
SUBARU observation (2003 January 26) does not allow the measurement of the faint object 
even if the measured total magnitude (V = 22.68 ± 0.015) suggests that we only see the 
ULX counterpart. The conclusion on the variability in the SUBARU observations is thus 
uncertain because of the photometric contamination by this fainter object. Nonetheless, all 
values are completely consistent with the luminosity range seen by the HST. We may be 
seeing similar variability as with the HST, but it is masked by the large uncertainties. 

An interesting characteristic of these HST observations is that they were coordinated 
with two simultaneous Chandra ones. Unfortunately, the observations are highly piled-up 
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preventing a detailed analysis. The derived fluxes suffer from a high uncertainty, with 
F x = 1-02+°;^ x KT 11 , LOOloJo x 10~ n erg s" 1 cm' 2 for the two observations. 
The optical variability observed AV ~ 0.14 mag corresponds to an increase of 14% in the 
optical flux. Considering the errors in the X-ray fluxes, it is obvious that nothing can be 
said about a possible correlation between the X-ray and optical fluxes. Detailed studies 
are needed to establish X-ray/optical correlations and it is clear that no conclusion can be 
drawn from only two data points. But the fact that X-ray and optical variations are seen 
in other ULXs which do not seem to be entirely related to the orbital movement of the 



companion star fjGrise et al. 



2008 . 



20091 ) are a strong motivation for future simultaneous 
X-ray/optical studies and will help to disentangle the mechanisms at the origin of the X-ray 
and optical emission. 



5. Conclusion 

We have used the high sensitivity of SUBARU/FOCAS and GEMIMI/GMOS-N 
and the high spatial resolution of HST/ ACS to make a detailed study of the optical 
counterpart of the ULX HoIX X-l and its immediate environment. Close to the ULX, 
we have highlighted a group of young stars spread out over 100 pc that appears to be an 
OB association rather than a gravitationally linked cluster. The photometric properties 
of these stars are consistent with the young population that surrounds them. This tight 
association is quite isolated, located at 1 kpc from the dense star-forming area of Holmberg 
IX. This recent and localized episode of star formation is possibly related to the past 
interaction between M81 and M82, similar to Holmberg IX and other structures formed at 
H I overdensities. We have estimated that the young stellar association linked to the ULX 
has an age < 20 Myr, a stellar-mass of ~ 1-3 xlO 3 M , and an upper mass limit for single 
stars of about 20 M & . We have newly estimated the reddening value near the association 
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and found a probable local extinction of E(B — V) ~ 0.18 mag in addition to the very 
modest Galactic reddening of E(B — V) ~ 0.08 mag. 

The optical counterpart of the ULX was identified by the presence of a broad, 
high-excitation He IIA4686 line that testifies to X-ray reprocessing in the system. The 
counterpart appears to be the most luminous object in the association. Using stellar 
evolutionary tracks, we have constrained the mass of the companion star to be > 10 M , 
even if the accretion disk contributes significantly to the optical luminosity. HoIX X-l 
displays a UV photometric excess compared to the surrounding population of stars. Its 
colors are not consistent with a star, and this prevents any spectral classification of the 
counterpart. The optical spectrum shows a blue continuum, consistent with the colors, 
with only high-excitation emission lines, similar to X-ray active LMXBs. Comparison of 
our photometric results with binary models suggest that the donor star may be as massive 
as ~ 30-50 M . 

Another interesting result of this work is the variability of the ULX counterpart, by 
~ 0.15 mag, clearly detected in the HST/ACS data set. However, we have insufficient 
information to determine if this variability comes from ellipsoidal variations, stochastically 
variable irradiation of the accretion disk or donor star, or some combination. 

Finally, although no stringent constraints can be put on the mass of the black hole from 
these data, we can, as for the ULX NGC 1313 X-2 (G08), emphasize that the formation of 
an intermediate-mass black hole (M > 100 M ) in the environment of HoIX X-l is very 
unlikely, given the stellar-mass of the association and its unbound nature. This is consistent 
with recent X-ray results showing that this source can be modeled by what appears to be 
super-Eddington emission, hence a small (< 100 M ) black hole. 
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54 53 52 

RIGHT ASCENSION 




54 53 52 

RIGHT ASCENSION 



Fig. 1. — : SUBARU/FOCAS composite images centered on the optical counterpart of the 
ULX. Left: image R (red) + V (green) + B (blue) ; right: image made of the nebular 
emission (H Q -continuum - red, [O IIl]-continuum - green) on which is superimposed the 
stellar emission coming from the V band (blue). 
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Table 1:: The SUBARU/FOCAS Observations for HoIX X-l 



Filter 


Exposure Time 

(s) 


Number of Exposures 


Observation Date 


Seeing 
(") 


B 


300 


2 


2003 Jan 25/26 


0.73/0.44 


V 


300 


4 


2003 Jan 25/26/26/Feb 09 


0.73/0.42/0.48/0.43 


R 


300 


2 


2003 Jan 25/26 


0.73/0.39 


H a 


900 


1 


2003 Jan 26 


0.51 


III 


900 


1 


2003 Jan 26 


0.53 


Spectrum 


1800/900 


2 


2003 Jan 26 


0.5 



18 19 20 21 22 23 24 25 

Subaru B 



18 19 20 21 22 23 24 25 

Subaru V 




18 19 20 21 22 23 24 25 

Subaru R 



Fig. 2. — : Photometric correlations between SUBARU data and the SDSS survey for three 
broadband filters (B, V, and R). The red line shows the range used to compute the photo- 
metric offset, using a weighted linear regression. 
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Table 2:: The HST/ ACS Observations for HoIX X-l 



ID 


Instrument 


Filter 


Date 


MJD at Mid-exposure 


Exposure Time (s) 


j8ola2010 


HRC 


F330W 


2004 Feb 7 


53042.51347 


2760 


j8ol02040 


WFC 


F435W 


2004 Feb 7 


53042.46162 


2520 


j8ol02030 


WFC 


F555W 


2004 Feb 7 


53042.42239 


1160 


j8ol02010 


WFC 


F8UW 


2004 Feb 7 


53042.41364 


1160 


j8o!06010 


WFC 


F555W 


2004 Mar 25 


53089.45604 


2400 
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Fig. 3.—: Difference between SUBARU/FOCAS and HST/ ACS photometry in the B (left 
panel) and V bands (right panel). For this comparison, we only used bright, isolated stars. 
We find that 50% of them lie within the combined photometric errors of the HST and 
SUBARU observations (solid red lines), and for most of the others the difference is only 
< 0.1 mag (dashed blue lines). 
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Table 3:: Luminosity of the ULX Counterpart in Different Filters from Our SUB- 
ARU/FOCAS Observations 

Filter Exposure Time (s) Date Magnitude Absolute Magnitude 

B 300 2003 Jan 26 22.72 ± 0.04 -6.1 

V 300 2003 Jan 26 22.71 ± 0.04 -5.9 

R 300 2003 Jan 26 22.64 ± 0.02 -5.8 



Table 4:: Brightness of the ULX Counterpart in Different Filters, from the HST/AGS Observations 



Filter 


Exposure Time 


Date 


Aperture Correction 


VEGAmag 


Johnson Magnitude 


F330W/U 


2760 s 


2004 Feb 7 


0.790 


20.833 ± 0.023 




FA35W/B 


2520 s 


2004 Feb 7 


0.427 


22.541 ± 0.015 


22.604 


F555W/V 


1160 s 


2004 Feb 7 


0.403 


22.632 ± 0.024 


22.609 


F8UW/I 


1160 s 


2004 Feb 7 


0.716 


22.334 ± 0.034 


22.328 


F555W/V 


2400 s 


2004 Mar 25 


0.384 


22.768 ± 0.013 


» 22.75 



Note. Magnitudes are expressed both in the HST/ ACS Vegamag system and in the Johnson- 
Cousins (UBVRT) system, when possible. 
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Table 5:: The Chandra Observations for HoIX X-l 



ID 


MJD 

Mid-exposure 


Exposure Time 
(ks) 


(10 21 cm" 2 ) 


r 


^^0.3-10 koV 

(10 -11 erg s _1 cm -2 ) 


Cash Stat./dof 


4751 
4752 


53042.45496 
53089.41938 


5 
5 


1.0 
1.0 


1 7 n+0.19 
1 74+0.19 

L - '^-0.19 


1 11+ 017 
i -- L1 -0.19 

i ir;+0.30 
i - ±o -0.21 


5.8/9 
6.6/7 



Table 6:: Position of the ULX, Derived from the Chandra and HST Data Sets and Compared 
with the XMM-Newton Position 



Telescope /Instrument 


Right Ascension 


Declination 


Chandra 


09 h 


57 m 


53! 32 


69°03'48'.T0 


HST/ ACS FA35W 


09 h 


57 m 


53!31 


69°03'48'/08 


HST/ ACS F555W 


09 h 


57 m 


53!31 


69 o 03'48'.'05 


HST/ ACS F8UW 


09 h 


57 m 


53!31 


69°03'48'.T0 


XMM-Newton 


09 h 


57 m 


53! 20 


69°03'47^7 



53.5 53.0 
RIGHT ASCENSION 



52.5 s 



53.5 53.0 
RIGHT ASCENSION 



Fig. 4.—: Identification of HoIX X-l ULX on an HST/ ACS image in the F435W filter 
(left) and in the F330W filter (right). The counterpart is designated by a blue cross (left). 
Slits from the SUBARU (position angle of 180°) and GEMINI observations (position angle 
of 90°) are overlaid on the right image. The Chandra (green circle) and XMM-Newton (red 
circle) positions are also overlaid with error circles of 0'/67 and (X'80 radius, respectively (90% 
confidence, including the error on the calibration of the optical image). We also note the 
association of stars to the east of the ULX. 



-33 - 



Table 7:: Observed and Dereddened Magnitudes and Colors of the ULX Coun- 
terpart, from the HST/ ACS Observations 



Filter /Magnitude 


Observed Magnitude/Color 


Dereddened Magnitude/Color 


B 


22.604 ± 0.015 


21.536 ± 0.015 


V 


22.609 ± 0.024 


21.780 ± 0.024 


I 


22.328 ± 0.034 


21.844 ± 0.034 


B-V 


-0.005 ± 0.028 


-0.247 ± 0.049 


V-I 


0.281 ± 0.042 


-0.064 ± 0.058 



Note. Values are expressed in the Johnson-Cousins (UBVRT) system. 




4600 4600 5000 5200 6000 6500 7000 



Wavelength (A) Wavelength (A) 

Fig. 5. — : One- dimensional FOCAS spectrum of the most luminous object in the field 
showing clearly the signature of the He II line at 4686 A, which proves that it is the optical 
counterpart of the ULX. The other annotated emission lines are residuals of the non-perfect 
subtraction process of the inhomogeneous nebula. 
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Continuum 

- - Hen 

— ► ULX 

— ► Other stars 



60 65 70 75 80 1140 1145 1150 1155 1160 1165 

Position (Pixel) Position (Pixel) 



Fig. 6. — : Spatial profiles of the continuum of the ULX (solid line) and of the He II emission 
line (dashed line) in the two co-added spectra (left: SUBARU data, right: GEMINI data). 
The continuum was taken between 4700 and 4714 A for the SUBARU data and averaged 
between 4662-4675 A and 4697-4711 A for the GEMINI data. The ULX is located near the 
peak at the pixel coordinate 71.5 (left) and near the pixel coordinate 1151.7 (right). The 
continuum of the ULX is slightly distorted, resulting from the presence of a nearby star at 
0'.'2 at the northwest of its position (on the right here). The vertical arrows show the position 
of the objects as derived from HST observations. In the SUBARU spectrum, we see that the 
peak of the He II emission line is coincident with the brightest object in the HST images. In 
the GEMINI spectrum, the position angle for these observations was set to 180° (90° for the 
SUBARU data). As can be seen in the HST image (Figure H]), three sources thus contribute 
to this spectrum. The tail of the emission at the left of the brightest object is due to a faint 
object located at / /6 south of the main object, the tail at the right is due to the closest 
object discussed already. Finally, we can see that the peak of the He II emission is again 
consistent with the position of the brightest object. This strongly suggests that the brightest 
object is the optical counterpart of the ULX. 
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4600 4700 4800 4900 5000 5100 



Wavelength (A) 

Fig. 7.—: 4600-5100 A part of the GEMINI/GMOS-N one-dimensional spectrum of HoIX X- 
1 optical counterpart, confirming the presence of the He II line at 4686 A . The other 
annotated line in emission is from the nebula in which the counterpart is located and which 
result probably from our rough subtraction considering the highly variable profile of the 
nebula. 
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Fig. 8. — : Location of the ULX compared with the denser star formation areas of Holm- 
berg IX, on the west edge, on an HST/ ACS observation in the F435W band. 
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Fig. 9. — : Left: true color image (blue: FA35W; green: F555W; red: F81AW) of the region 
around HoIX X-l, from HST/AGS observations. 1" represents 17.4 pc at the distance of 
Holmberg IX. The stellar environment is quite poor in the region close to the ULX. We can 
mainly see a blue association localized to the east of the counterpart. Right: same view of 
the vicinity of the counterpart, in the F435W filter. The ULX counterpart is the brightest 
point-like source at the center of the image. Contours of the H a emission (at 10%, 20%, 
40%, 60% and 80% flux level above the background) are overplotted, from our SUBARU 
observations. The stellar association is localized or projected inside the H a nebula. 
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Fig. 10. — : Spectral energy distribution (SED) for HoIX X-l and six of the following brightest 
stars that have been shifted down for clarity purpose. Upper panels: the use of the Galactic 
extinction E(B — V) = 0.08 fails to reproduce the slope of the expected OB stars. Lower 
panels: increasing the extinction to E(B — V) = 0.26 makes the SEDs of the stars fairly 
consistent with those of OB stars templates. Left panels show only the SEDs of the stars 
near the ULX (triangles) and right panels show in addition stellar templates designated by 
their spectral type. 



-40- 
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Fig. 11. — : HST/ ACS color- magnitude diagrams for the stellar field (2000 x 2000 pixels, i.e., 
100" x 100") around the ULX. HST/ACS magnitudes were transformed into the Johnson- 
Cousins system. Padova isochrones for stars of different ages are overplotted. Typical 
photometric errors are also plotted. Data have been corrected for the derived extinction 
{E(B — V) = 0.26 mag), the bar at the top left corner illustrating this effect. "Stars 
near ULX" are stars not located inside the H a nebula but within 25" of the ULX position. 
"Candidate cluster stars" are those located inside the nebula and detected in the UV image 
(F330W filter). "Probable cluster members" are those inside the nebula but not detected 
in the UV image, with B — V < 0.5. The "Other UV stars" are stars detected in the UV 
but which are located outside the nebula. Finally, "Field stars" are objects which are not 
compatible with any of these criteria. Left panel: color-magnitude diagram in the (B,V) 
system. Right panel: color- magnitude diagram in the (V,I) system. The same isochrones 
are plotted in the two panels, i.e., for 1, 5, 10, 20, 50, 200, and 500 Myr at Z = 0.008. 
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Fig. 12. — : HST/ ACS color-magnitude diagrams (Johnson- Cousins photometry) with 
Padova evolutionary tracks for stars of different initial masses (4, 7, 10, 15, 20, 30, and 
40 M Q ) with Z = 0.008. The legend is the same as Figure [TTJ Left panel: color-magnitude 
diagram in the (B, V) system. Right panel: color-magnitude diagram in the (V, I) system. 
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Fig. 13. — : HST/ ACS color-magnitude diagram for the stellar field around the ULX in the 
synthetic (U,B) system. The ULX counterpart (red diamond) is clearly bluer than the other 
stars of the association, contrary to the other color-magnitude diagrams. 
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Fig. 14. — : Spectral energy distribution (SED) for two ULXs including the object of our 
study, with the SEDs of different OB star templates. The two dot dashed black lines are 
power laws with indices of 2.0 and 1/3, respectively corresponding to a Rayleigh- Jeans tail 
and the expected spectrum of a thin accretion disk. 
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Table 8:: Magnitudes of the ULX Counterpart and of its Close Object in the 
V Filter Derived from SUBARU/FOCAS and HST/ ACS Observations 



Telescope 


Date 




UuLX 


^close object 


Vtotal 


SUBARU/FOCAS 


2003 Jan 


25 










22.62 


± 


0.02 


SUBARU/FOCAS 


2003 Jan 


26 


22.71 ± 0.037 


24.81 


± 


0.26 


22.56 


± 


0.26 


SUBARU/FOCAS 


2003 Jan 


26 


22.68 ± 0.015 














SUBARU/FOCAS 


2003 Feb 


9 


22.67 ± 0.025 


24.49 


± 


0.12 


22.48 


± 


0.12 


HST/ACS 


2004 Feb 


4 


22.609 ± 0.024 


24.54 


± 


0.03 


22.44 


± 


0.04 


HST/ACS 


2004 Mar 


25 


22.745 ± 0.013 


24.50 


± 


0.02 


22.55 


± 


0.02 



Note. Magnitudes are expressed in the Johnson-Cousins system (UBVI). 
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